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Abstract 


The Norwegian Lundehund is a highly endangered native dog breed. Low fertility and high 
frequency predisposition to intestinal disorder imply inbreeding depression. We assessed 
the genetic diversity of the Lundehund population from pedigree data and evaluated the 
potential of optimal contribution selection and cross-breeding in the long-term management 
of the Lundehund population. The current Norwegian Lundehund population is highly inbred 
and has lost 38.8% of the genetic diversity in the base population. Effective population size 
estimates varied between 13 and 82 depending on the method used. Optimal contribution 
selection alone facilitates no improvement in the current situation in the Lundehund due to 
the extremely high relatedness of the whole population. Addition of (replacement with) 10 
breeding candidates of foreign breed to 30 Lundehund breeders reduced the parental addi- 
tive genetic relationship by 40-42% (48-53%). Immediate actions are needed to increase 
the genetic diversity in the current Lundehund population. The only option to secure the con- 
servation of this rare breed is to introduce individuals from foreign breeds as breeding 
candidates. 


Introduction 


The Norwegian Lundehund is a highly endangered dog breed, native to Norway [1], listed as a 
national Norwegian dog breed. The name of the Norwegian Lundehund originates from its 
function to retrieve puffins alive (Fratercula arctica, lunde in Norwegian) from their nests in 
screes and burrows on steep mountainsides along the Northern Norwegian coast [1,2]. There 
is evidence that this hunting method dates back as far as to the 16" century [1]. Lundehund 
were highly appreciated and secured both income (down feather) and food for families living 
in fisherman’s villages near seabird colonies [1]. In the 1850s hunting with nets became popu- 
lar and the need for Lundehund disappeared. Consequently, the number of Lundehund 
reduced drastically. At the outbreak of the Second World War the population size was 50 indi- 
viduals [1]. Since then the Lundehund population experienced two severe bottlenecks. Asa 
result the current Lundehund population stems from only a very few related individuals. 
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It is expected that the bottlenecks have caused a great loss of genetic diversity in the Lunde- 
hund population. Visual evidence for this loss of diversity is the reduced variation of coat col- 
ours. Before the Second World War black and white Lundehund existed, whereas at present 
these variants are extinct [1]. Reduced litter size has been suggested as a sign of inbreeding 
depression in the Lundehund. Currently, the Norwegian Lundehund club reports average litter 
sizes of 2.8. Low average litter size is a reflection of a high proportion of one-puppy litters and 
near complete absence of litters with six or more puppies [3]. Moreover, breeders have 
reported problems relative to fertility, such as “invisible” heat, behavioural problems when 
mating, both in males and females, and low sperm quality. Overall, reduced litter size most 
likely reflects inbreeding depression in several traits, both physiological and behavioural, as 
well as a high frequency of non-viable embryos due to lethal alleles. 

In the 1960s breeders discovered that many dogs were suffering from intestinal problems, 
later diagnosed as intestinal lymphangiectasia (IL) [1]. Untreated IL is a life threatening dis- 
ease, but even with an early diagnosis and veterinary treatment the survivors risk a relapse. 
About 30% of the reported causes of death in the Lundehund in Norway are due to IL [3]. Due 
to the complicated disease mechanisms and unknown mode of inheritance [4], breeding 
efforts have failed to reduce the high frequency of this disease in Lundehund [1]. It may be 
speculated that such a high disease frequency is also a consequence of inbreeding. Pfahler and 
Distl [5] suggested that selection for body size or male fertility might have predisposed the 
Lundehund to gastroenteropathies. 

Molecular genetic studies have verified extremely low genetic diversity in the Lundehund. 
Melis et al. [2] reported observed heterozygosity of 0.075 for 26 microsatellite loci in a Norwe- 
gian sample of Lundehund. Based on SNP markers (single nucleotide polymorphism), Pfahler 
and Distl [6] found an observed heterozygosity of 0.047 and high Fj; (0.87) indicating close 
relatedness between the individuals. In comparison, the heterozygosity in five Danish dog 
breeds was estimated to be between 0.27 and 0.36 [7]. Similarly, a recent study verified signifi- 
cantly lower genetic diversity in the Lundehund than in three other Nordic Spitz breeds: Nor- 
wegian Buhund, Icelandic sheepdog and Norrbottenspets (unpublished data, Stronen A). 

The Norwegian Lundehund club, the organization responsible for breeding and conserva- 
tion of the Lundehund, has initiated a cross-breeding project to improve health and welfare of 
the population. In particular, it is critical to reduce the risk of developing IL and improve fertil- 
ity without compromising the preservation of the distinctive characteristics of the breed. The 
Norwegian Kennel Club, the head organization for breed organizations, has granted permis- 
sion to cross three phenotypically similar breeds into the Lundehund: Norwegian Buhund, Ice- 
landic sheepdog and Norrbottenspets. To date, only a few attempts of crossbreeding have been 
made, and only two of those matings have resulted in litters. 

The sustainable management of small populations is dependent on the long term control of 
increases in additive genetic relationships. Optimal contribution selection (OCS) is a method 
that allows maximisation of genetic merit of a cohort of animals, while setting constraints on the 
average relationships between the individuals. OCS has proven to be a robust method also in 
complex breeding scenarios, imposed by various restrictions, realising most of the long-term 
genetic gain obtained by OCS without restrictions [8]. Many small native populations have a 
pure conservation focus, and consequently the genetic contributions of selection candidates, 
individuals available for breeding, are optimised only based on the additive genetic relationships. 

Our overall goal was to evaluate the conservation possibilities of the current Lundehund 
population. Consequently, we assessed the genetic diversity of the Lundehund population 
from pedigree data in terms of levels of inbreeding and relationships, probability of gene ori- 
gin, and effective population size. We furthermore evaluated the potential of OCS and cross- 
breeding in the long-term management of the Lundehund population. 


PLOS ONE | DOI:10.1371/journal.pone.0170039 January 20, 2017 2/16 


®-PLOS | ONE 


Cross-Breeding to Conserve the Norwegian Lundehund 


Materials and Methods 
Data 


Pedigree data from 1930 to 2015 comprising records of 5433 individuals was obtained from 
the Norwegian Lundehund Club. The data included identification of the individual (number 
and name), identification number of sire and dam, country of origin, country of residence, 
sex, birth date, information whether the dog is alive or not, and number of offspring. A total of 
1220 individuals were registered as alive, of which 50.4% were living in Norway. The largest 
populations of Lundehund outside Norway are situated in Sweden (N = 156), Finland 
(N = 106) and Denmark (N = 78). The distribution of the current global Lundehund popula- 
tion is presented in Table 1. For OCS analyses dogs born between 2007 and 2014 were consid- 
ered as breeding candidates. To avoid that parent(s) and offspring were allocated to the same 
reproductive cohort the contemporary group (later referred to as a time step) was defined as a 
six-month period within each birth year, e.g. 2014/1 (January-June) and 2014/2 (July- 
December). 

The current pedigree dates back to the 1930s and these founder dogs were considered as 
unrelated. If the dogs used to construct the current population were related, the level of relat- 
edness and average inbreeding in the current population will be underestimated. 


Statistical analyses 


Population analyses relative to inbreeding, additive genetic relationships and pedigree com- 
pleteness were carried out using EVA 2.1 [8,9]. The effective population sizes were estimated 
by ENDOG v.4.8 [10]. ENDOG utilizes, by default, either the full (F) or restricted (R) pedi- 
gree data set, where individuals with unknown parentage and generation coefficient < 2.0 
were discarded from a restricted pedigree data set, in order to take into account that 


Table 1. Distribution of the current global Lundehund population. 


Country Number of dogs % 

Norway | 617 50.41 
Sweden | 156 12.75 
Finland | 106 8.66 
Denmark 78 6.37 
The Netherlands 67 5.47 
Germany 55 4.49 
USA 49 4.00 
Switzerland 30 2.45 
France 18 1.47 
Czech Republic 17 1.39 
Austria 8 0.65 
Luxembourg 7 0.57 
Belgium 3 0.25 
Iceland 2 0.16 
Cyprus 2 0.16 
Poland 2 0.16 
Croatia 1 0.08 
Lithuania 1 0.08 
Spain 1 0.08 
Total 1220 100 


doi:10.1371/journal.pone.0170039.t001 
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Table 2. Summary of the methods used for estimation of effective population size. 


Method Definition of AF Regression on | Pedigree data set 
la b Number of complete generations |F 
ie) b Maximum number of generations F 
Ic b Complete generation equivalents | F 
ll CED Complete generation equivalents R 
Ixb ] 
ll Fey Birth years F 
IVa =) Birth years | F 
IVb a2) Complete generation equivalents R 
V AF, =1—-'\/1—F, Not relevant IR 
Vi AG = 1? 87 = 6, Not relevant 'R 


F = full pedigree data set; R = restricted pedigree data set with individuals with unknown parentage and generation coefficient (complete generation 
equivalent) < 2.0 discarded; V and VI definition is for individual increase in inbreeding and co-ancestry, respectively. Full description of methods III, IV b, V 
and VI are found in [11], [12], [13,14] and [15], respectively. 


doi:10.1371/journal.pone.0170039.t002 


ancestral information has to be known for a minimum of two generations back to detect 
inbreeding. Individuals with a low level of pedigree information would contribute to overes- 
timation of N.. ue effective population size based on different definitions of AF was com- 
puted as N, = =; 


F* 
First, three simple regression based methods were used for estimating AF from the full data 


set (Table 2). Individual inbreeding coefficients were regressed on the number of complete 
generations (method I a), the maximum number of generations (method I b) and the number 
of equivalent generations traced back for each individual (method I c) and the regression coef- 
ficient (b) was used as the increase in inbreeding between two generations. Consequently, N 
was estimated as N, = =; [10]. The regression coefficient (b) of the ee inbreeding coef- 


ficients over the equivalent generations was used to approximate AF ~ in method II. In 


= ae =) 
the third regression based method (method III), following Gutiérrez et al. [11], the inbreeding 
coefficients are regressed on birth years and the increase in inbreeding between two genera- 
tions is defined as the proate of regression coefficient (b) and the average generation interval 
(D: lx b. From this, AF = ra 
question. Additionally, values of log(1 — F;) were regressed on birth years (method IV a) or 
(1-e) 
i 


where F; is the average inbreeding of the population in 


complete generation equivalents (method IV b) (Perez-Enciso 1995) and AF = was 
calculated. 
In method V, N, via individual increases in inbreeding (AF;) was estimated according to 


Gutiérrez et al.[13]. “Realized effective size” was computed as N, = sap Where AF is computed 


— i 

as AF = ( by =) / , where F; is the individual inbreeding coefficient, t the generation 
N 

coefficient of individual i, and N is number of individuals in the reference population [13,14]. 


The standard error of N, can be computed as: ox- = 2N,”ox¢ ro [14]. Finally, a method 


based on the increase in co-ancestry for all pairs of individuals in the reference population was 
used for the estimation of N, (method VI). Averaging out the increase in co-ancestry in all 


pairs of individuals in the reference population provides an estimate of Ac and N,. = =: The 


standard error of the estimate N,. is calculated in a similar manner as in the method of Gutiér- 
rez et al. [13,14]. See, Cervantes et al. [15] for a complete description of this method. 
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Parameters derived from the probabilities of gene origin were used as complementary to 
increase in inbreeding and effective population size to describe genetic diversity in the Lunde- 
hund population. The effective number of founders or the founder equivalent (f,) is defined as 
the number of equally contributing founders that would be expected to result in the same 
genetic diversity as in the population under study, thus accounting for the loss in genetic vari- 
ability due to unequal founder contributions [16]. It is calculated as an inverse of the sum of 
the squared proportional contributions (p,) of the founders f, = '/ 3 p> The effective num- 


ber of ancestors (f,) is defined as the minimum number of ancestors (not necessarily founders) 
explaining the complete genetic diversity of the population. Parameter f, is calculated as f, but 
replacing p; with marginal ancestral contributions (m;), and thus unlike f,, accounts for poten- 
tial bottlenecks in the pedigree [17]. Finally, the founder genome equivalent f,. [16,18] takes 
into account the number of founder alleles already lost from the population both due to bottle- 
necks and genetic drift. This parameter accounts for all causes of gene loss and describes how 
many alleles and in which frequencies they have been maintained in a given locus [17]. It may 
be approximated using both the proportion of genes contributed by a given founder in the 
descendant population and the expected proportion of alleles retained within the descendant 
population [18]. Lacy [18] introduced the exact definition of f,, which may be obtained from 
the average co-ancestry of a predefined reference population tas f,, = '/,, [19]. The effective 


-l1 
number of non-founders was calculated as f,, = [ = 4] , measuring the amount of genetic 
fe Se 


drift since the foundation of the population. The parameters f,, f,- and f,,. were used to derive 
measures of genetic diversity. The degree of genetic diversity relative to the base accounting 
for loss of diversity due to unequal founder contribution was calculated as GD* = 1 — ra [18] 


and accounting for unequal founder contribution, bottlenecks and genetic drift as GD = 
1— oe = 1— 4% [18,19]. The difference GD* — GD = a. measures the loss in diversity due to 
“Tge “Ine 


drift accumulated over non-founder generations. Using the above expressions, the propor- 


tional loss of diversity due to random genetic drift was simplified to 1 — =1- a Conse- 


quently, proportional loss due to unequal founder contribution is ‘e, Reference population in 
the estimation of parameters based on the probabilities of gene origin was defined as all living 
individuals (Table 1). 

The software package ENDOG was used to obtain the N, estimates, but all, except that of 
Cervantes et al. [15], are easily calculated from the general output files of both EVA and 
ENDOG following the formulas given above and the full descriptions of the methods in the lit- 
erature [11,13-15]. Parameters f, and f, were obtained from ENDOG whereas f,. was calcu- 
lated as explained above from EVA output files. 

EVA 2.1 software was used for the OCS analyses [8,9]. To assess the possible effect of 
importation on the management of the Norwegian Lundehund population, three different 
data sets were analysed: global, Nordic and Norway (Table 3). Additionally, following the 


Table 3. Number of male and female candidates and the relatedness between male and female candidates in the different OCS scenarios. 


Males Females | Relatedness* 
Global | 462 478 0.77 
Nordic 365 | 362 | 0.77 
Norway 237 | 224 | 0.77 
Norway, no-crossbreds 286 | 220 | 0.78 


* average genetic relationship between male and female candidates including relationship to itself (1+F;). 


doi:10.1371/journal.pone.0170039.t003 


PLOS ONE | DOI:10.1371/journal.pone.0170039 January 20, 2017 5/16 


oe 
(6) . PLOS | ONE Cross-Breeding to Conserve the Norwegian Lundehund 
oe 


1.0 

0.9 

0.8 

0.7 ue rl 

O in a ee 

De) ----PCI7 
EA ——PCI10 
] 

0.5433 
id 
i‘ 

0.4 +4 

0.3 

0.2 

1964/1 1969/1 1974/1 1979/1 1984/1 1989/1 1994/1 1999/1 2004/1 2009/1 2014/1 

Time step 


Fig 1. Pedigree completeness (PCI) between 1964-2015 for pedigree completeness indices of 5 (PCI5), 7 (PCI7) and 10 (PCI10) 
generations. Time Steps are defined as 6-month periods within a year: /1 = January-June, /2 = July-December. 


doi:10.1371/journal.pone.0170039.g001 


recommendations of the Norwegian Lundehund club prohibiting the use of cross-bred individ- 
uals in breeding at the moment, an additional analysis excluding cross-bred individuals was car- 
ried out. Three scenarios of limited use of males within each data set were tested, allowing males 
to be used in 1, 5 or 10 matings. A total of 20 matings were optimized in the OCS analyses. 


Results and Discussion 
Pedigree completeness 


Pedigree completeness [20] of the Lundehund from the start of 1964 is presented in Fig 1. The 
most conservative pedigree completeness index (PCI), looking back 10 generations (PCI10), 
reached 90% completeness during the early 1990s (Fig 1). The same level of PCI was reached 
in the early 1980s and late 1970s for PCI7 and PCI5, respectively. During the same time period, 
the depth of the pedigree increased from 6 to 19 generations. 

The drops in PCI in 2004 and 2014 are due to the proportion of unknown parentage: 9% 
(2004) and 12% (2014). Two incidences of drops in the average level of inbreeding coincide 
with these declines in PCI10. The first case is a result of “unknown” parentage of an imported 
individual causing an underestimation of the level of inbreeding in that cohort, as all Lunde- 
hund originate from Norway and thus are related to each other. The second case, in 2014, is a 
real fall in inbreeding as a consequence of using a Norwegian Buhund dam in breeding. High 
PCI of the Lundehund pedigree data ensures reliable estimation of the average levels of 
inbreeding in the current population. 


PLOS ONE | DOI:10.1371/journal.pone.0170039 January 20, 2017 6/16 


®-PLOS | ONE 


Cross-Breeding to Conserve the Norwegian Lundehund 


Genetic contribution of founders and ancestors 


There were 49 founder animals (f) in the data set. If an animal with one unknown parent was 
defined as a half founder, the actual base population consisted of 43.5 individuals. Out of 
these, only 11 contributed to the last cohort of 2015 (percentage of genetic contribution in 
parentheses). Eight of the founders (7.0%) are parents born in the 1930s, each having one con- 
tributing offspring. One male founder (4.7%) has three offspring born in the 1960s, of which 
two have further contributed to the last cohort. The founder male contributing most to the last 
cohort has 10 offspring all born in the 1960s, which all have further contributed to the popula- 
tion. This individual is the second most contributing ancestor in the pedigree, with 35% con- 
tribution to the last cohort. 

The ancestor with the highest contribution in the pedigree is a female with 18 offspring 
born in the 1960s. Her contribution to the last cohort is 41%. Her pedigree can be traced back 
to the base population through three rounds of full-sib mating, making her 50% inbred. The 
cumulative proportion of genetic variation explained by these two most contributing ancestors 
is 76%. In other words, on average 76% of alleles in individuals born in the first half of 2015 
originate from these two ancestors. 


Relatedness and levels and trends of inbreeding 


The average co-ancestry and inbreeding coefficient in the last fully registered cohort 
(2014/2) was 35.6% and 33.9%, respectively. The development of the level of inbreeding (F) 
between 1964/1 and 2015/1 together with PCI10 is shown in Fig 2, panel A. The estimated 
rate of inbreeding (AF) throughout the pedigree is 1.0% per generation if linear regression 
of F over time was used for estimation. AF based on the individual increases in inbreeding 
was 4% [13]. 

In total, 33 individuals registered in the pedigree had an inbreeding coefficient equal to or 
larger than 0.5. One generation of self-fertilisation or three generations of full-sib mating result 
in inbreeding coefficients of 0.5. Of the total number of matings, 0.98% (N = 53), 2.04% 

(N = 111) and 0.81% (N = 44) were full-sib, half-sib or parent-offspring matings, respectively. 

The global Lundehund population is highly inbred. To our knowledge, this is the second 
highest level of inbreeding reported for dogs in the scientific literature. Since 1996 the aver- 
age within cohort inbreeding has been relatively stable between 36-38%, with the exception 
of two cohorts (2004/1; 2014/2), where importation and the use of a foreign breed lowered 
the average inbreeding to 34% (Fig 2, panel A). In the Lundehund, the expected and observed 
levels of inbreeding have not significantly deviated since the early 1970s (alpha in in Fig 2, 
panel B). This implies that the intended avoidance of mating between close relatives have 
had no effect on the levels of inbreeding as the co-ancestry between all breeding candidates is 
uniform. 

Several studies have reported average inbreeding coefficients for various dog breeds 
either assessed by genealogical or molecular data. Only one pedigree-based estimate of 
inbreeding for Lundehund has been reported prior to our study; Pfahler and Distl [5] 
have recently estimated inbreeding in Lundehund from pedigree data reaching back 11 
generations. Their estimate of mean inbreeding of 10% is likely an underestimate as there 
was no knowledge on relatedness of the founders of the pedigree. Typically, low or moder- 
ate pedigree-based average inbreeding coefficients has been published in the scientific liter- 
ature for numerous dog breeds (a selection summarized in Table 4). Only a few studies 
have reported levels of average inbreeding comparable to those found in this study (0.21- 
0.45). In some cases, level of inbreeding is likely to be underestimated due to a shallow or 
sparse pedigree. 
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Fig 2. Average inbreeding and pedigree completeness in the Lundehund population between time steps 1964/1 
and 2015/1. (A) Average inbreeding (F) and pedigree completeness considering 10 generations (PCI). (B) Expected and 


observed inbreeding and their deviation (alpha). Time steps are defined as 6-month periods within a year: /1 = January- 
June, /2 = July-December. 


doi:10.1371/journal.pone.0170039.g002 
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Table 4. Average inbreeding coefficients in various dog breeds, a selection. 


Breed Mean F Reference 
Rough Collie 0.073 [21] 
Labrador Retriever 0.024 [21] 
Greyhound 0.058 [21] 
Polish Hound 0.37 [22] 
“Healthy” breeds 0.045* [23] 
“Unhealthy” breeds 0.025* [23] 
Norrbottenspets 0.047 [23] 
Norwegian Buhund 0.054 [23] 
The Pyrenean Mountain Dog 0.04 [24] 
Barbet 0.124 [24] 
The Pyrenean Shepherd 0.088 [25] 
The Czechoslovakian Wolfdog 0.002 [25] 
Great Dane 0.044 [25] 
Nova Scotia Duck Tolling Retreiver 0.26 [26] 
Lancashire Heeler 0.10 [26] 
Icelandic Sheepdog 0.21 [27] 
Bouvier des Ardennes 0.45 [28] 


*dogs alive at the end of 2010, mean of healthy (N = 11) or unhealthy breeds (N = 16) 


doi:10.1371/journal.pone.0170039.t004 


Effective population size 


Table 5 summarizes the results of the estimation of N, with different methods. Regression 
based methods estimated up to over 6-fold N, (34-82) compared to estimates based on indi- 
vidual increase in inbreeding or co-ancestry in the Lundehund (13) [13,14]. As long as the full 
pedigree information is utilised, methods accounting for full pedigree history were relatively 
invariant from the choice of the sample used for N, estimation, whereas the regression on 
birth date was very sensitive to the choice of the time period (Table 6). When regression based 
estimates were obtained using a time period with relatively stable level of inbreeding, overesti- 


mation of N, was drastic (Fig 2, Table 6). 


Table 5. Estimates of effective population size (N,) estimated with different methods. For methods V 


and VI, standard error in parentheses. 


Method Ne 
Whole pedigree 
la: | Regression on complete generations 34 
|b: | Regression on maximum generations 82 
Ic: | Regression on equivalent generations 47 
Ill: | Regression on birth date 33 
IV a: | Log regression on birth date 37 
Restricted pedigree 
Il: | Regression on equivalent generations 38 
IV b: | Log regression on equivalent generations 41 
V: | Individual increase in inbreeding 13 (1.35) 
VI: | Individual increase in co-ancestry 13 (0.45) 


doi:10.1371/journal.pone.0170039.t005 
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Table 6. Effect of the choice of reference population on the estimation of effective population size 
(Ne). 


Ne 
Time period Individual increase in inbreeding Regression on bir date 
1964-2015 13 47 
1930-1972 8 13 
1973-1999 | 11 34 
2000-2015 | 16 | 134 


Individual increase in inbreeding refers to the method of Gutiérrez et a/. [12,13] and regression on birth date 
to that of Gutiérrez et al. [11]. 


doi:10.1371/journal.pone.0170039.t006 


N, limits are often used as indicators for sustainability and endangerment. Populations 
with lower N, than 100 are said likely to be at risk for substantial losses in genetic diversity, 
and those with an N, lower than 50 at high risk of detrimental effects of inbreeding [29]. No 
earlier pedigree based estimates of N. are available for Lundehund. The present study estimates 
the N, of the global Lundehund population to be very low (13). Based on the number of breed- 
ing animals, the effective population size for the Lundehund has earlier been reported to be 
approximately 80 [3]. Pfahler and Distl [5] estimated inbreeding and N, from a molecular data 
of 28 individuals of Lundehund. Their estimate for the last 7-200 generations was 10-13, in 
good agreement with our pedigree based estimate, supporting a genetically small and closed 
population [5]. 

Regression based and highly variable estimates of N, have been reported for a large group 
of dog breeds [21,23-26,29]. Leroy et al. [24] reported realized effective population sizes of 20- 
147 in nine French dog breeds, and in a more recent paper an N, of 46-2136 in 61 dog breeds 
in France [25]. According to Calboli et al. [21] N, varied between 17 (Greyhound) and 114 
(Labrador retriever) in10 dog breeds registered in the UK Kennel Club since 1970. In a more 
recent paper, Lewis et al. [29] reported trends in genetic diversity in all dog breeds currently 
recognized by the UK Kennel Club. Most likely the estimates therein are overestimates, as the 
time span of this research is 1980-2014 and relatedness before 1980 was not accounted for. A 
recent study by Wijnrocx et al. [28] reported highly variable N, estimates based on both indi- 
vidual increase in inbreeding (3.2-829.1) and individual increase in co-ancestry (4.5-334.9) in 
total of 23 Belgian dog breeds. The variability in the reported estimates of N, arises from the 
breed history, the amount and quality of the pedigree information, and in particular, from the 
method used to obtain the estimates. In regression based methods, inbreeding coefficients are 
regressed on either generations or year of birth. Consequently, the choice of the time period 
used for calculation of N, becomes decisive: flat or negative trends of inbreeding will result in 
illogical estimates of N, [14]. As long as the full pedigree information is utilised methods based 
on individual increases in inbreeding/co-ancestry describe the breed history, inclusive effects 
such as mating policy, drift and selection. Thus, these methods are expected to result in esti- 
mates invariant from the choice of the sample used for N, calculation [14]. Additionally, this 
approach allows ascertainment of the confidence in the estimate, because the corresponding 
standard error of the estimate is easily computed. 


Probability of gene origin 


Parameters of probability of gene origin were estimated to be 6, 3, 1.3 and 1.6 for f., fa» foe and 
Jne- Only a small proportion of the base population contributed to the current Lundehund pop- 
ulation, reflected by low ratio of f, to f (0.12). Comparison of f, and f, in Lundehund indicates 
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substantial loss of genetic variation due to bottlenecks. The parameter f, is not very useful for 
assessing genetic diversity. First, it only reflects the inequalities in the proportional contribu- 
tions of founders and does not account for genetic drift [16], and second, it becomes constant 
after a few generations [19] as the genetic contributions of founders will converge. The param- 
eter f,. is diminished by both unequal proportional contribution of founders and by genetic 
drift, that is, discounts f, to account for the irreversible losses of founder alleles [16]. The value 
of fge represents the cumulative loss of genetic diversity since the base population and directly 
relates to N, [18,19]. Ratios of the above mentioned parameters of gene origin have been used 
to interpret the loss in genetic diversity by several authors [23,24,30], but definition and inter- 
pretation of the parameter ratios have been diverse and non-quantifiable (in absolute sense). 

Genetic diversity measures GD, GD* and GD*-GD were used to quantify the impacts of 
unequal founder contribution, bottlenecks and random genetic drift. It was estimated that 
8.3% of the variation in the base population was lost due to unequal founder contribution 
(GD* = 0.917) and additional 30.5% due to genetic drift accumulated over non-founder popu- 
lations (GD*-GD). Consequently, 38.8% of the genetic diversity in the base population was 
lost from the current living Lundehund population. Proportionally, 78.5% of the total loss was 
due to random genetic drift and 21.5% due to unequal founder contributions. In fact, the total 
loss in genetic diversity is simply the average co-ancestry of the chosen reference population 
[19], which has been between 0.35 and 0.40 since 1980s in Lundehund. During that same time, 
PCI5 has improved from 0.94 to 1.0, except for the time steps including imported or crossbred 
individuals. Consequently the f,. has been constant, as has the total loss in genetic diversity. 

Molecular genetic studies have showed extremely low genetic diversity in the Lundehund 
[2,6,31] and close relatedness between individuals [6]. In a recent study (unpublished data, 
Stronen A) SNP genotypes were used to estimate heterozygosity in the Lundehund population. 
The level of heterozygosity was significantly lower (0.038) than in other Nordic Spitz breeds 
(0.230-0.298). Using this value and our estimate of the level of inbreeding, the heterozygosity 
prior to the bottlenecks is estimated to be 0.061. This, again, indicates that the Lundehund was 
already a genetically small population at the time of the start of the pedigree recordings, and 
that the founder dogs were likely highly related. 


Optimal contribution selection 


Very similar results of OCS were obtained from the data sets global, Nordic and Norway, and 
only the latter is presented here (Table 7). Relaxing the limitation of number of matings 
allowed per male results in favouring only males least related to the female candidates. In the 
case of 10 matings allowed per male, only two cross-bred males born in 2014 were selected as 
sires. If only one mating per male was allowed, the average relatedness of the selected males to 


Table 7. Optimal contribution selection from the Norwegian Lundehund data, 20 matings requested. 


Maximum number of matings per male 


a 1 | 5 10 
Average relationship _ [ov [oa | 
Average inbreeding 0.34 | 0.29 0.21 
Generation interval 10.0 8.3 6.2. 
AF -0.11 | -0.22 -0.35 
Number of males 20 | 7 2 
Number of females 20 | 20 20 
Maximum avoidance of inbreeding 0.31 0.24 0.21 


doi:10.1371/journal.pone.0170039.t007 
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the female candidates was 72%, whereas the relatedness between all male candidates to all 
female candidates was 77%. When 5 or 10 matings per male were allowed, the average related- 
ness of the selected males to the female candidates was 64% and 38%, respectively. The reduced 
generation interval with more matings allowed per male, is due to the more intense use of the 
young crossbred males. 

When the individuals of cross-bred origin (total of 6 individuals) were excluded as breeding 
candidates, OCS had no effect. The different scenarios of restricted male use resulted in identi- 
cal levels of average relatedness and inbreeding in the next generation as observed in the last 
cohort, 76% and 38%, respectively. 

The results of the population analysis of Lundehund quantify the narrow genetic base well 
documented in the breed’s history [1]. Although dog breeders have been very aware of avoid- 
ing the mating of close relatives this has had very little effect in the Lundehund population due 
to the extremely high relatedness among all individuals, which is a result of a low genetic diver- 
sity already early in the population history, combined with several severe bottlenecks. This is 
illustrated by identical levels of expected (under random mating) and observed inbreeding 
over the past decades. 

Due to the extremely high relatedness, OCS alone facilitates no improvement in the current 
situation in the Lundehund. This was clearly demonstrated by invariant OCS solutions over 
regions (global, Nordic, Norway). If cross-bred individuals were excluded as breeding candi- 
dates, no favourable changes in relatedness were achieved independent on the number of mat- 
ings allowed per male. 

In situations where genetic variation is depleted in a population, migration is the only 
method to introduce new variation into the population and enable sustainable management 
and selection [29]. Based on the results of this study, together with previously published 
results on Lundehund genetic diversity [2,6,31], it is clear that immediate actions are needed 
to increase the genetic diversity in the current Lundehund population, and that the only 
option to secure the conservation of this rare breed is to introduce individuals from foreign 
breeds. 

The initialisation of cross-breeding does not offer an immediate remedy to the serious 
problem of relatedness and inbreeding in the Lundehund population. How much the additive 
genetic relationship of the breeding candidates, half of the expected F of the progeny, is 
decreased by the introduction of foreign breeds depends on the number of breeders, number 
of breeds and individuals of foreign breeds introduced, within-breed additive genetic relation- 
ships, and the average inbreeding coefficients of all breeds. For Lundehund (30 purebred 
breeding candidates with an average additive relationship, a, of 0.77), introduction of addi- 
tional 10 parents of foreign breed reduced the average parental a by approximately 40%, down 
to 0.42-0.48 (Fig 3). This was achieved independently of whether one, two or three breeds 
(with equal proportions, within breed a 0.08-0.4, and F = a/2) were introduced. If instead of 
adding new candidates, ten purebred individuals were replaced by parents from other breeds, 
a reduction of 48-53% in parental a was observed. When breeders of foreign breeds were 
added to the default Lundehund parental group, the number of breeds and different levels of a 
and F of the foreign populations seemed to have a minor effect on the parental a, whereas the 
number of individuals introduced was decisive (Fig 3). When Lundehund were replaced by 
foreign breed individuals, larger differences were observed between scenarios of different 
number of breeds, F and a (Fig 3). The effect of cross-breeding on the next generation’s relat- 
edness is conditional on the proportional contributions of foreign breed and purebred parents, 
which in turn depends on the breeding decisions of private dog owners. The conservation of 
Lundehund allele variants is difficult to predict as it is conditional on the future use of the 
cross-bred individuals. 
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Fig 3. Effect of introduction of foreign breeds into the default Lundehund parental population (N = 30, additive 
genetic relationship 0.77). Default parental population: N = 30, additive genetic relationship (a) 0.77, inbreeding a/2. (A) 
Introduction of one breed. (B) Introduction of two or three breeds. 


doi:10.1371/journal.pone.0170039.g003 
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Conclusions 


The Lundehund is an endangered Norwegian dog breed with high recreational and cultural- 
historical value. Current population is showing signs of inbreeding depression, due to 
extremely high relatedness of the whole population. Breeding optimisation, including OCS, 
offer no remedy unless individuals of foreign breed are introduced to the Lundehund. Our 
data strongly suggests that crossbred individuals should immediately be accepted as part of the 
main population and as breeding candidates. 
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